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Abstract 



We review progress in calculating one-loop scattering ampli- 
tudes required for next-to-leading-order corrections to QCD pro- 
cesses. The underlying technical developments include the spinor 
helicity formalism, color decompositions, supersymmetry, string 
theory, factorization and unitarity. We provide explicit examples 
illustrating these techniques. 



1.1 Importance of Diagrammatic Calculations 

Gauge theories form the backbone of the Standard Model. The weak- 
coupling perturbative expansion of gauge theory scattering amplitudes, 
carried out by means of Feynman diagrams, has led to theoretical predic- 
tions in remarkable agreement with high-energy collider data [Q. This 
high-precision agreement places strong bounds on new physics. In the 
strong-interaction sector of the Standard Model — described by quan- 
tum chromodynamics — the precision is not as great as in the elec- 
troweak sector. QCD is asymptotically free, so the strong coupling 
constant as becomes weak at large momentum transfers, justifying a 
perturbative expansion Physical quantities do depend on nonper- 
turbative, long-distance QCD, in the form of quantities such as parton 
distribution and fragmentation functions, as well as on the physics of 
hadronization. In many processes at modern colliders, however, the 
dominant theoretical uncertainties are due to an incomplete knowledge 
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of the perturbation series, rather than to our relative ignorance of non- 
perturbative aspects of scattering processes. The situation is exacer- 
bated by the slow approach to asymptopia {as is of order 0.1 at the 
100 GeV scale), and by the presence of large logarithms of ratios of 
scales. 

The leading-order (LO) term in the as expansion of a QCD cross- 
section conies simply from squaring a tree-level scattering amplitude. 
Efficient techniques for computing QCD tree amplitudes have been avail- 
able for some time now Q , and the results have provided a basic theoret- 
ical description of QCD processes and thereby estimates of QCD back- 
grounds to new physics searches. Unfortunately, higher order correc- 
tions, especially those enhanced by logarithms, can be sizeable. The ul- 
traviolet logarithms manifest themselves in the residual renormalization- 
scale dependence of a finite-order prediction. The renormalization scale 
fj-R is introduced in order to define the coupling constant; renormaliza- 
tion group invariance requires any physical quantity to be independent 
of it. However, when a perturbative expansion is truncated at a finite 
order, residual /ifl-dependence appears, because the cancellation takes 
place across different orders in as ■ Calculations at next-to-leading order 
(NLO) in as significantly reduce the dependence on fin as compared to 
leading order. As an example, fig. |l] shows the comparison of the LO 
and NLO theoretical predictions to the experimental measurement of a 
point in the single-jet inclusive distribution. Note the good agreement 
between NLO theory and experiment and the significant reduction of 
theoretical uncertainties, compared to the LO calculation. 

Infrared logarithms arise because jet processes involve more than 
one scale, at the very least a scale characterizing the jet size in addition 
to the hard scale of the short-distance scattering, and because of the 
infrared divergences of perturbative QCD. These divergences transform 
the perturbation expansion for such quantities from one in as alone 
to one in aslog^ ym and as log ym in addition to as, where ym is a 
jet 'resolution' parameter. All three must be small for the perturbation 
expansion to be reliable; but the first two cannot be calculated in an LO 
calculation. Only in an NLO calculation are the corresponding terms 
determined quantitatively, and only at this order can one establish the 
reliability of the perturbative calculation. 

Beyond the logarithmically-enhanced corrections, the 0{as) correc- 
tions to most jet observables are larger than non-perturbative power 
corrections and corrections due to quark masses, and are thus the most 
important ones to calculate in order to refine the precision of theoretical 
predictions. 
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Figure 1: The inclusive cross section for single-jet production in pp 
collisions at ^/s = 1.8 TeV and jet transverse energy Et = 100 GeV 
(using MRSDg structure functions showing the sensitivity of the 
LO result to the choice of renormalization scale, fiR, and the reduced 
sensitivity at NLO. The CDF data shown is extracted from ref. the 
band shows statistical errors only. 



Despite the need for higher-order QCD computations, at present 
no quantities have been computed beyond next-to-next-to-leading or- 
der (NNLO), and the only quantities that have been computed fully 
at NNLO are totally inclusive quantities such as the total cross-section 
for e~^e~ annihilation into hadrons, and the QCD corrections to various 
sum rules in deeply inelastic scattering 0. At NLO, there are many 
complete calculations (in the form of computer programs producing nu- 
merical results) for a variety of processes, but at present results are 
still limited to where the basic process has four external legs (counting 
electroweak vector bosons rather than their decay products as external 
legs) . The following are examples of calculations which are relevant for 
current experiments but have not yet been performed or assembled: 

1. NLO corrections to three-jet production at hadron colliders. These 
contributions would allow a measurement of as (via the three-jet 
to two-jet ratio) at the highest experimentally available momen- 
tum transfers, as well as next-to- leading-order studies of jet struc- 
ture. 

2. NLO corrections to -I- multi-jet production at hadron coUid- 



4 



ers. These processes form a background to the t quark signal at 
Fermilab. 

3. NLO corrections to e+e^ 4 jets. At the Z resonance, this is the 
lowest-order process in which the quark and gluon color charges 
can be measured independently. It will also be useful for ruling 
out the presence of light colored fermions (or scalars). At LEP2 
it is a background to threshold production of W pairs, when both 
W's decay hadronically. 

4. NNLO corrections to e'^e~ 3 jets. These corrections are the 
dominant uncertainty in a precision extraction of from hadronic 
event shapes at the Z j^. 

In any of these processes, deviations of experimental results from the 
theoretical predictions could indicate new physics. 

Why do these higher-order QCD corrections remain uncalculated? 
NLO corrections can be divided into real and virtual parts. (See fig. |^.) 
Real corrections arise from the emission of one additional parton into 
the final state, and are straightforward to compute from tree amplitudes 
with one more leg than the LO tree amplitude. Virtual corrections arise 
from the interference of the LO tree amplitude with a one-loop ampli- 
tude. Each contribution is infrared divergent, but the divergences can- 
cel in the sum, after integrating the real contribution over "unobserved" 
partons in the final state and factorizing initial state singularities 
into the definition of parton distributions in an incoming hadron [l0[| . 
The remaining finite integrations are typically performed with a numer- 
ical program ||Tl|] . 

While the numerical evaluation of NLO corrections can be non- 
trivial, the major analytical bottleneck is simply the availability of one- 
loop amplitudes, which enter into the virtual corrections. In partic- 
ular, one-loop amplitudes with more than four external legs (and all 
quarks massless), which are required for the higher-order corrections 
listed above, have only recently become available, thanks to the devel- 
opment of new calculational techniques. The purpose of this review is 
to provide an introduction to some of these techniques, together with 
worked-out examples. 

Our emphasis will be on obtaining compact analytic results. In 
general, it is preferable to have such results for matrix elements, even 
though they are ultimately inserted into numerical programs for com- 
puting cross-sections. Without compact results, numerical instabilities 
can arise from the vanishing of spurious denominators in the expression. 



5 



2 



2 



\ 



/ 



/ — \ 



Figure 2: In (a) the parton subprocesses required for the LO contribu- 
tion to two-jet production at hadron coUiders are shown schematically. 
In (b) the corresponding real and virtual NLO contributions are shown. 

With analytic forms it is also easier to compare independent calcula- 
tions, to understand better how to organize calculations, and even to 
obtain results for an arbitrary number of external legs |l3|, |l^ |l^ . 

1.2 Difficulty of Brute- Force Calculations 

Gauge theories have an elegant construction based on the principle of 
local gauge invariance. The QCD Lagrangian for massless quarks q is 



where the covariant derivative = — igA^/\/2 and field strength 
F^j/ = iV2[-D/j, -Di/]/<7 are given in terms of the matrix-valued gauge 
connection = AJ^T°.0 Since Cqcd depends on a single couphng 
constant g, all the interactions are dictated by gauge symmetry. Un- 
fortunately, the Feynman diagram expansion does not respect this in- 
variance, because the quantization procedure fixes the gauge symmetry. 
Individual diagrams are not gauge invariant, and are often more com- 
plicated than the final sum over diagrams. The non-abelian gluon self- 
interactions coming from the cubic and quartic terms in eq. (Q) have a 
complicated index structure and momentum-dependence. So while it is 
straightforward in principle to compute both tree and loop amplitudes 
by drawing all Feynman diagrams and evaluating them, in practice this 
method becomes extremely inefficient and cumbersome as the number 
of external legs grows. For five or more external legs there are a large 

^ The normalization Tr(r''T'') = (5"* of the fundamental-representation genera- 
tors accounts for the \/2's here; it serves to eliminate the \/2's from the partial 
amplitudes defined below. 
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Figure 3: The five-gluon pentagon diagram. 



number of kinematic variables, which aUow the construction of comph- 
cated expressions. Indeed, intermediate expressions tend to be vastly 
more complicated than the final results, when the latter are represented 
in an appropriate way. 

As an example consider the five-gluon pentagon diagram, depicted 
in fig. |[ which would be encountered in a brute-force computation of 
NLO corrections to three-jet production at a hadron collider. Each of 
the five non-abelian three-point vertices in the diagram is given by 

V^^'^{k,p, q) = r''{vupip - q)^ + VpM - k). + VpAk - P)p) , (2) 

where /"'"^ are the SU{3) structure constants, fc, p and q the momenta, 
and rjfii, the Minkowski metric. As the non-abelian vertex has six terms, 
a rough estimate of the number of terms is about 6^ . Each term is asso- 
ciated with a loop integral which evaluates to an expression on the order 
of a page in length. This means that one is faced with about 10'' pages 
of algebra for this single diagram. As bad as this brute-force approach 
might seem, the situation is actually worse, because of the structure of 
the results. After evaluating the integrals and summing over a few hun- 
dred more diagrams one obtains expressions of the form where 
the factors Ni are polynomials in the gluon polarization vectors and ex- 
ternal momenta, and the Di (polynomials in the external invariants) are 
produced when the loop integrals are reduced to a standard set of func- 
tions. In general the Di contain spurious kinematic singularities which 
cancel only after combining many terms over a common denominator; 
this causes an explosion of terms in the numerator. 

In contrast to the complexity of intermediate expressions, the final 
results can be strikingly simple. F or e xample, the five gluon amplitudes 



which we shall describe in section 3^ are remarkably compact. 
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1.3 Non-traditional Approaches 



Substantial progress has been made in the past decade in improving 
the calculation of tree- level amplitudes. Four ideas which have played 
an important role are the spinor helicity method for gluon polarization 
vectors |16 , the color decomposition supersymmetry identities 

, and the Berends and Giele recurrence relations . Although these 
ideas form a basis for the one-loop techniques described here, they have 
been extensively reviewed in ref. ^, |2l| , and we permit ourselves only a 
brief review below. 

As illustrated by the pentagon example above, one-loop computa- 
tions are significantly more complicated than tree computations, so fur- 
ther techniques are useful for preventing an explosion in algebra. The 
additional ideas which we shall discuss in this review involve string the- 
ory, supersymmetry, unitarity, and factorization. String theory, for ex- 
ample, suggests better gauge choices, a supersymmetric decomposition 
of amplitudes, and an improved disentanglement of color and kinemat- 
ics. Approaches based on unitarity and factorization make use of the an- 
alytic properties of amplitudes to build further amplitudes using known 
ones. Since these approaches use gauge-invariant quantities as the basic 
building blocks of new amplitudes they tend to be extremely efficient. 
Although we shall not discuss recursion relations here, Mahlon has made 



considerable progress in applying these to one- loop amplitudes 1 13 . For 



simplicity, we demonstrate the methods for amplitudes where all exter- 
nal particles are giuons, even though most of the techniques (or analogs 
of them) can be applied to amplitudes with external fermions as well. 

To date, these techniques have allowed for the computation of all 
one-loop five-parton helicity amplitudes |22 , ^|, a s well several in- 



finite sequences of one-loop amplitudes |l2| , 13| , l4 15| . The five-parton 
amplitudes are currently being incorporated into numerical programs for 
NLO three-jet production at hadron colliders, the first item on the list 
in Section 1.1 j25|. Thus the analytical bottleneck to NLO corrections 
is yielding to the new techniques described in this review. 



2 PRIMITIVE AMPLITUDES 

In this section, we briefiy review the use of color and helicity information 
to decompose amplitudes into 'primitive amplitudes'. These building 
blocks have a much simpler analytic structure than the full amplitudes, 
a fact which will be exploited in subsequent sections. We also review 
the application of supersymmetric Ward identities to QCD. 
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2.1 Color Decomposition 

Color decompositions have a long history, dating back to Chan-Paton 
factors in early formulations of string theory |l7| . They are also related 
to the "double-line" formalism introduced by 't Hooft in the large-iVc 
(number of colors) approach to QCD ||2^, although here we will not 
make any large- A^c or "leading-color" approximations. The basic idea is 
to use group theory to break up an amplitude into gauge-invariant pieces 
which are composed of Feynman diagrams with a fixed cyclic ordering of 
external legs. These pieces are simpler because poles and cuts can only 
appear in kinematic invariants made out of cyclicly adjacent sums of 
momenta, of the form (fc^ -f /ci+i -f • • • -|- /cj)^. At the four-point level this 
is not so important, because only one of the three Mandelstam variables 
s,t,u is thereby excluded; but as the number of external legs grows, the 
total number of invariants grows much faster than the number of cyclicly 
adjacent ones. The following brief review focuses on results needed later, 
rather than derivations. A more complete discussion can be found in 
refs. 0,|,|27[0|4[|l|. 

We first generalize the gauge group of QCD to SU{Nc), with the 
quarks transforming in the fundamental representation. The simplest 
way to implement the color decomposition in field theory is by rewriting 
the group structure constants appearing in Feynman diagrams in terms 
of fundamental representation matrices 

After making this substitution in a generic Feynman diagram we obtain 
a large number of traces, many sharing y's with contracted indices, of 
the form Tr(. . .T" . . ) Tr(. . . . . .) ... Tr(. . .). If external quarks are 

present, then in addition to the traces there will be some strings of T"'s 
terminated by fundamental indices. To reduce the number of traces and 
strings to a minimum, we rearrange the contracted T^'s, using 

^^)n^a^n ^ s^5^~^5i^5^, (4) 

where the sum over a is implicit. (If all lines are in the adjoint rep- 
resentation the second term drops out by a U{1) decoupling identity 
1^, |l^, which follows from the lack of a 'photon' self-coupling.) A par- 
tial amplitude is the coefficient of a given color trace in the resulting 
color decomposition of the amplitude. 
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For example, in the n-gluon tree amplitude, application of eq. (Q) 
reduces all color factors to single traces. Thus its decomposition is 

A'r = g''-' E Tr(a(l)...f7(n))ArKl),...,^H), (5) 

Tes„/z„ 

where A^''"'^ are the partial amplitudes, Tr(l ...n) = Tr(r''i • • 
with ai the color index of the i-th external gluon, and S'„/Z„ is the set 
of non-cyclic permutations of {1, 2, . . . , n}, corresponding to the set of 
inequivalent traces. The labels on the gluon momenta ki and polariza- 
tion vectors e^, implicit in eq. (j^), are also to be permuted by a. In the 
next subsection we will go over to a helicity basis, and the label i will 
be replaced by i^' , with Ai the (outgoing) gluon helicity. Similarly, tree 
amplitudes with a pair of external quarks can be reduced to a sum over 
single strings of matrices, (T"^ • • •T°")j^\ and so on. For a proof that 
individual partial amplitudes are gauge invariant, see ref. 

At one loop, additional color structures are possible; in the n-gluon 
amplitude double traces appear as well as single traces. For example, 
the color decomposition of the one-loop five-gluon amplitude is 

E A^cTrKl) . . . a(5))A5;i(a(l), . . . , (t(5)) 
+ E TrKlM2))Tr(a(3)a(4)a(5))A5;3(a(l),a(2);a(3),a(4),a(5)) ; 

0-eS5/(5'2XS3) 

(6) 

as in eq. the permutation sums are over all inequivalent traces. For 
gauge group U{Nc), the partial amplitudes A5-2 multiplying traces of 
the form Tr(l)Tr(2345) would also have to be included, but for SU{Nc) 
the trace of a single generator vanishes. The decomposition of the n- 
gluon amplitude into single-trace (An-i) and double-trace {An-j>2) com- 
ponents is entirely analogous. Were one to consider the large- iVc limit, 
the single-trace terms would give rise to the leading contributions, and 
we will refer to the corresponding partial amplitudes as leading-color 
partial amplitudes; the double-trace terms have subleading-color partial 
amplitudes as coefficients. 

The rules for constructing leading-color partial amplitudes such as 
All°° and An-i are color-ordered Feynman rules, which are depicted in 
fig. ^ for the standard Lorentz-Feynman gauge. These rules are ob- 
tained from ordinary Feynman rules by restricting attention to a given 
ordering of color matrices. Applying eq. (j^) to eq. (|^) and extract- 
ing the coefficient of t^y{T°'T^T'') gives the color-ordered three-vertex 



.l-loop _ 5 2e 
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Figure 4: Color-ordered Feynman rules in Lorentz-Feynman gauge. 
Curly lines represent gluons and lines with arrows fermions. 



in fig. ^; and similarly for the color-ordered four-vertex, the coefficient 
of ^y{T''T^T''T'^). The only diagrams to be computed are those that 
can be drawn in a planar fashion with the external legs following the 
ordering of the color trace under consideration. 

The immediate advantage of rewriting Feynman rules in this way 
is that fewer diagrams contribute to a given partial amplitude, and its 
analytic structure is simpler. As a simple example, with conventional 
Feynman diagrams one would have a total of four conventional Feyn- 
man diagrams, depicted in fig. ^for the four-point tree amplitude. With 
color-ordered Feynman rules one would compute the partial amplitude 
yl4(l,2,3,4) associated with the color trace Tr(r°iT"=T''3T"^), omit- 
ting diagram ||c since the ordering of the legs do not follow the ordering 
of the color trace. Thus ^4(1,2,3,4) has no pole in (fci -I- ^3)^. The 
other partial amplitudes can be obtained by permuting the arguments 
of ^4(1, 2, 3, 4). For the five-gluon amplitude, there are 10 color-ordered 
diagrams as opposed to 40 total. Obviously the simplifications obtained 
using partial amplitudes increase rapidly with the number of external 
legs. 

At one loop, one also has to compute subleading-color partial ampli- 
tudes, such as the double-trace coefficients ^5;3 in eq. (^, which cannot 
be obtained directly from color-ordered rules. Fortunately there exist 
general formulas relating such quantities to permutation sums of color- 
ordered objects ^ For example, the gluon-loop contribution 
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Figure 5: The four-point Feynman diagrams. Color-ordered Feynman 
rules do not include diagram (c) for ^4(1, 2, 3, 4). 

to the four-gluon amplitude can be found from the relation 

A4;3(l, 2; 3, 4) = A4;i(l, 2, 3, 4) + Ai,i{l, 3, 2, 4) + ^4;i(2, 1, 3, 4) 

+ A4;i(2,3,l,4) + A4;i(3,l,2,4) + A4;i(3,2,l,4). ^ 

Such formulae can be derived from string theory (28[ although the 
most straightforward way to prove them is using color flow diagrams in 
field theory . To understand formula ^ heuristically, it is useful to 
focus on the box diagram. Using ordinary Feynman rules and expanding 
out the structure constants using eqs. (||) and it is straightforward 
to check that the box diagrams contribute to A^-^i and A^-^^ in such a 
way that eq. (0) is satisfied. Roughly speaking, gauge invariance then 
requires the remaining diagrams to tag along properly with the box 
diagram. 

Thus we can restrict our discussion henceforth to amplitudes with 
a fixed ordering of external legs, which we call primitive amplitudes. 
In the n-gluon cases discussed above, the set of primitive amplitudes 
coincides with the leading-color partial amplitudes A^^'^^ and but 
this is not always the case. For example, one-loop amplitudes with 
external fermions have leading-color (as well as subleading-color) partial 
amplitudes that are sums of several primitive amplitudes [M . 



2.2 Spinor Helicity Formalism 

In explicit calculations, it is very convenient to adopt a helicity (circu- 
lar polarization) basis for external gluons. The spinor helicity formal- 
ism expresses the positive- and negative-helicity polarization vectors 
in terms of massless Weyl spinors |fc^), 
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where q is an arbitrary null 'reference' momentum which drops out of the 
final gauge-invariant amplitudes. (Changing q is equivalent to perform- 
ing a gauge transformation on the external legs.) We use the compact 
notation 

(fcr|fc+) ^(zj), {kt\k-)^M. (9) 

These spinor products are crossing-symmetric, antisymmetric in their 
arguments, and satisfy 

{ij)[j{\^2h-kj=s,j. (10) 

Helicity amplitudes can be given a manifestly crossing symmetric rep- 
resentation, with the convention that a helicity label corresponds to an 
outgoing particle; the helicity of an incoming particle is reversed. As 
we shall discuss in Section ^ in the coUinear limit where ki and kj be- 
come parallel, helicity amplitudes have a square-root singular behavior, 
^ --^= ^ -jj-jj ~ whose magnitude and phase are captured con- 
cisely by the spinor products. This helps explain why spinor products 
provide an extremely compact representation of amplitudes. 
In performing calculations, the Schouten identity is useful, 

{zj){kl) = {tl){kj) + {tk){jl) . (11) 

A more complete discussion, including further identities and numerical 
representations of the spinor products, can be found in refs. ||l^, ||, [2l| |. 

To maximize the benefit of the spinor helicity formalism for loop 
amplitudes we must choose a compatible regularization scheme. In 
conventional dimensional regularization p9[ |, the polarization vectors 
are (4 — 2e)-dimensional; this is incompatible with the spinor helicity 
method, which assumes four-dimensional polarizations. To avoid this 
problem, we modify the regularization scheme so all helicity states are 
four-dimensional and only the loop momentum is continued to (4 — 2e) 
dimensions. This is the four-dimensional- helicity (FDH) scheme ||2^ , 
which has been shown to be equivalent ||30|| to an appropriate helicity 
formulation of Sicgel's dimensional-reduction scheme |3l| at one-loop. 
The conversion between schemes has been given in ref. ]30| , so there is 
no loss of generality in choosing the FDH scheme. 

2.3 Parity and Charge Conjugation 

The reader might worry that the color and helicity decompositions will 
lead to a huge proliferation in the number of primitive or partial ampli- 
tudes that have to be computed. In fact, this does not happen, thanks 
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to the group theory relations mentioned above, plus the discrete symme- 
tries of parity and charge conjugation. Parity simultaneously reverses 
all helicities in an amplitude; eq. (^) shows that it is implemented by 
the exchange (ij) ^ [j i]. Charge conjugation is related to the anti- 
symmetry of the color-ordered rules; for pure-glue partial amplitudes it 
takes the form of a reflection identity, 

4r(l,2,...,n) = (-1)" 2,1). (12) 

For amplitudes with external quarks, it allows one to exchange a quark 
and anti-quark, or equivalently to flip the helicity on a quark line. 

As an example, with the use of parity and cyclic (Z5) symmetry, we 
can reduce the five-gluon amplitude at tree level to a combination of 
just four independent partial amplitudes: 

^■--(l+,2+,3+,4+,5+), 4™(l-,2+,3+,4+,5+), 
4-'=(l-,2-,3+,4+,5+), <°<=(l-,2+,3-,4+,5+). 

Furthermore, the first two partial amplitudes here vanish (see below), 
and there is a group theory (C/(l) decoupling) relation between the last 
two 1^, so there is only one independent nonvanishing object to cal- 
culate. At one loop there are four independent objects — eq. (|lj) with 
^trcc replaced by A^-i — but only the last two contribute to the NLO 
cross-section, due to the tree-level vanishings. The explicit expression 



for A5;i(l ,2 , 3+, 4+, 5+) is given in section 3.2 



2.4 Supersymmetry Identities 

What does supersymmetry have to do with a non-supersymmetric the- 
ory such as QCD? The answer is that tree-level QCD is "effectively" 
supersymmetric | p^ . Consider an n-gluon tree amplitude. It has no 
loops in it, so it has no fermion loops in it. The fermions in the the- 
ory might as well be in the adjoint representation, that is, the theory 
might as well be a super Yang-Mills theory. Pure-gluon tree amplitudes 
in QCD are indeed identical to those in the supersymmetric theory, 
and are thus related by supersymmetry to amplitudes with fermions 
(the gluinos). It is however more useful to think of such relations as 
connecting partial amplitudes. These relations arc the so-called super- 
symmetric Ward identities (SWI) |2l|. They connect pure-gluon 
partial amplitudes to partial amplitudes with a quark pair, because after 
the color information has been stripped off, the latter are identical to 
partial amplitudes with gluinos instead of quarks. Using the SWI saves 
computational labor [p^. 
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The SWI relate amplitudes with all external gluons, g, to amplitudes 
where a pair of gluons is replaced by a pair of gluinos, A, or a pair of 
complex scalars, 0. Specifically, the SWI that we shall make use of in 
later sections are 

ASUSY(a-,5+ ...,5+_^,A+)=0, 
^f'('/'r,52+,---,5+-i,0+)=O, 

^SUSY/ » - + - »+N _ O''^) aSVSY, - + - +^ 

V^l I 52 : ■ ■ • J i/j T ■ ■ ■ i-'^n ) — JJl^ " Iffl I 52 ' ■ • • ' yj 1 ■ ■ • ' 5n ^ 1 

(14) 

where '. . .' denotes positive-helicity gluons, and the helicity assignments 
on (f) refer to particle or antiparticle assignments rather than genuine he- 
licity. At tree level, these identities hold for all QCD partial amplitudes; 
for supersymmetric partial amplitudes, they hold to all orders in per- 
turbation theory. 

At one loop QCD "knows" that it is not supersymmetric, but one can 
still perform a supersymmetric decomposition of a QCD amplitude (see 



section 3.2), for which the supersymmetric components of the ampli- 
tude will obey eq. (|lj). One may also use the identities to find relations 
amongst non-supersymmetric contributions. For example, in iV = 1 
super- Yang-Mills, one can use the first of the identities in eq. (|lj) to de- 
duce that fermion- and gluon-loop contributions are equal and opposite 
for n-gluon amplitudes with maximal helicity violation. By considering 
an = 2 theory, with one gluon, two gluinos and one (complex) scalar, 
one deduces that all three types of loop contribution must be propor- 
tional to each other. We therefore obtain, for SU{Nc) QCD with Ug 
massless complex scalars and nf massless Dirac fermions, 

A„;i(5^g+, . . . ,3+) = (l + ^ - ^) A-f '-(fff ,52+, • • ■ ,5+) , (15) 

where ^1^':^''"' is the contribution of a single scalar and the factors of 
1/Nc are the conversion factors between the adjoint and fundamental 
representation loops. (Note that adjoint representation complex scalars 
have two states, but we have chosen the normalization that scalars in 
the fundamental representation — Nc(BNc — have four states, the same 
as for their would-be superpartner fermions.) 
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Figure 6: A single string diagram implicitly contains all field theory 
Feynman diagrams. 

3 STRING-INSPIRED METHODS 

String theory has provided a number of improvements in the calcu- 
lation of one-loop amplitudes. Originally, we used it to derive a set 
of diagrammatic computational rules for calculating gluon amplitudes 
, |3^. Such rules were used in the first computation of the one-loop 
five-gluon amplitudes |Q. The same methods also work well for grav- 
ity calculations |Q . One of the authors has reviewed the string-based 
rules in ref. [Q and we shall not do so here. Other approaches to the 
string-based rules have been formulated In particular, the first- 

quantized particle world-line has led to a rather efficient computation of 
the two-loop QED /3-function |^ and of coefhcients of high-dimension 
operators in effective actions. 

3.1 String Organization 

The basic motivation for the use of string theory follows from the com- 
pact representation it provides for amplitudes: at each loop order there 
is only a single closed string diagram. As depicted in fig. ^, the string 
theory diagram contains within it all the Feynman diagrams, including 
contributions of the entire tower of superheavy string excitations. The 
unwanted superheavy contributions are removed by taking the "low- 
energy limit" where all external momentum invariants are much less 
than the string tension. This limit picks out different regions of in- 
tegration in the string diagram (see fig. |^), corresponding roughly to 
particle-like diagrams, but with different, string-based, rules js^. 

Given knowledge of the string-based rules and organization, one may 
also formulate a conventional field-theory framework which mimics them 
1^ (at least for one- loop multiparton amplitudes), but which can be ap- 
plied more broadly (for example, to amplitudes with external fermions). 
At one loop, key ingredients of this string-inspired framework are: use 
of a special gauge which is a hybrid of Gervais-Neveu gauge |Q and 
background- field gauge |Q; improved color decompositions; systematic 
organization of the algebra; and a second-order formalism for fermions 
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Figure 7: The color-ordered Gervais-Neveu gauge three- and four-point 
vertices. 




1| which helps make supersymmetry relations manifest. 



Gervais-Neveu gauge, originally derived from the low-energy limit 
of tree-level string amplitudes , has the following gauge- fixed action 
(ignoring ghosts): 

5°^ = I d^x(--T,\Fl\-\^,\{^■A-^gA^|^2f-^. (16) 

The color-ordered Feynman rules derived from this action are depicted 
in fig. 1^; comparing them to the color-ordered vertices for the standard 
Lorentz-Feynman gauge (fig. we see that the three-point and four- 
point vertices have, respectively, half and a third as many terms, showing 
why the Gervais-Neveu gauge is simpler for tree-level calculations. 

Given this understanding of the string reorganization of tree am- 
plitudes, one might guess that string theory would best be described 
by the Gervais-Neveu gauge at one loop as well. However, the gauge 
most closely resembling the string organization of one-loop amplitudes 
is a hybrid gauge involving both background-field and Gervais-Neveu 
gauges Q . To quantize in a background-field gauge Q one splits the 
gauge field into a classical background field and a fluctuating quantum 
field, = A^ + A^ , and imposes the gauge condition D^A'j^ = 0, 
where = — -^gA^ is the background-field covariant derivative, 

with A^ evaluated in the adjoint representation. The Feynman-gauge 
version of the gauge-fixed action is (again ignoring ghosts), 

^B-^sd^ |d^:,(^-lTT[F'^J~^TT[i^■AQ~^g[A^,AQ]/V2f]). (17) 

The color-ordered background-field gauge vertices which arise from ex- 
panding eq. dr^ ) are depicted in fig. |^. Here we show only the vertices 
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Figure 8: The color-ordered back] 
and four-point vertices. Dashed hne 
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;round-field Feynman gauge three- 
i represent either ghosts or scalars. 



bihncar in the quantum field A^. These suffice for computing the one- 
loop effective action r[A^], since A'^ describes the gluon propagating 
around the loop while A^ describes a gluon emerging from the loop. 

Any one-loop diagram can be split into a one-particle-irreducible 
(IPI) part, or loop part, along with a set of tree diagrams sewn onto 
the loop. Now, r[A-^] is invariant with respect to A^ gauge transfor- 
mations Q. Therefore we may use any single gauge to compute the 
trees which are to be sewn onto the IPI parts of the diagrams. Indeed, 
the string-motivated recipe is to use background-field gauge only for 
the IPI or loop vertices, and Gervais-Neveu gauge for the remaining 
tree vertices |^^. This approach retains the above- noted advantages of 
Gervais-Neveu gauge for tree computations, while avoiding the compli- 
cated ghost interactions this nonlinear gauge would entail if it were used 
inside the loop. The advantage of the background-field gauge inside the 
loop is that the loop momentum appears in only the first term in the 
tri-linear gauge vertex in fig. ^ the last two terms contain only the ex- 
ternal momentum k. (In general, the most complicated loop integrals to 
evaluate are those with the most insertions of the loop momentum in the 
numerator.) Furthermore, the first term matches the scalar-scalar-gluon 
vertex, up to the ri,jp factor. Thus in background-field gauge the lead- 
ing loop-momentum behavior of one-particle-irreducible graphs with a 
gluon in the loop is very similar to that of graphs with a scalar in the 
loop. Note also that the interactions of a scalar and of a ghost with the 
background field are identical, up to the overall minus sign for a ghost 
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loop. In the next subseetion we elaborate further on these relations. 



3.2 Supersymmetric Decomposition 

String theory suggests a natural decomposition of QCD amplitudes into 
supersymmetric and non-supersymmetric parts. For example, for an n- 
gluon one-loop amplitude the contributions of a fermion and of a gluon 
circulating in the loop can be decomposed as 



^fermion yjscalar i aN—1 

^n:l — ^n:l > ^n:l 



/■gluon _ /(scalar a aN=1 , a 
^n;l — ^n;l ~ *^n;l + ^n;l • 

Here the "scalar" superscript denotes the contribution of a complex 
scalar in the loop; the = 1 superscript refers to the contribution 
of a iV = 1 supersymmetric chiral multiplet, consisting of a complex 
scalar and a Weyl fermion; and the N = 4 label refers to a vector 
supermultiplet, consisting of three complex scalars, four Weyl fermions 
and a single gluon, all in the adjoint representation. (We have assumed 
the use of a supersymmetry preserving regulator ||3^, HJ, ^ in these 
equations, and the vector-multiplet loop is defined to include the ghost 
loop.) 

The two supersymmetric components of eq. ( p^ ) have important can- 
cellations in their leading loop-momentum behavior. The simplest way 
to see this is via the scalar, fermion and gluon loop contributions to the 
background-field effective action, 

^'^^'"[A] =lndetjo[ {D^) , 

ptcrmionj^] = ^ hi detf//^, {D^ ~ g^CJ^^'^ ^) , (IQ) 

pgiuonj^j = Indctf,]'/' _ gSM-_p^^/V2) + lndet[o] {D'') , 

where ^u^i/ and Yj^y are respectively the spin-i and spin-1 Lorentz 
generators, and where det[j] is the one-loop determinant for a particle 
of spin J in the loop. The fermionic contribution has been rewritten in 
second-order form using 



Indet^/^ 



1% W-5lndetf/;,, (#) , 



(20) 



In an m-point IPI graph, the leading behavior of each contribution 
in eq. (M) for large loop momentum I is The leading term always 
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comes from the term in eq. (|T^) , because F^^ contains only the exter- 
nal momenta, not the loop momentum. Using Tr[o](l) — 1, Tr[]^/2](1) — 
Tr[i](l) = 4, we see that the term cancels between the scalar and 
fermion loop, and between the fermion and gluon loop; hence it cancels 
in any supersymmetric linear combination. Subleading terms in super- 
symmetric combinations come from using one or more factors of F in 
generating a graph; each F costs one power of £. Terms with a lone F 
vanish, thanks to Tr cr^j, = TrE^^ — 0. This reduces the leading power 
in an m-point IPI graph from i™ down to This argument can be 

extended to any amplitude in a supersymmetric gauge theory ||l5| and is 
related to the improved ultraviolet behavior of supersymmetric ampli- 
tudes. For the amplitude ^^f*, a comparison of the traces of products 
of two and three cr^,y's and E^^'s shows that further cancellations reduce 
the leading power behavior all the way down to £"^~*. This result can 
also be derived by superspace techniques In a gauge other than 

Feynman background-field gauge, the cancellations involving the gluon 
loop would no longer happen diagram by diagram. 

We illustrate the supersymmetric decomposition with the five-gluon 
amplitude, ^5.1(1^, 2~, 3+, 4+, S"*"), whose components (|T8|) are p2| 



\N=4 



crA' 



4^ 



+ In f —fi^) In f + 



1 _ 1 [ln(-S23) +ln(-S5i)] +2 



icT (1 2)^ ((2 3) [3 4] (4 1) + (2 4) [4 5] (51 

S51 - S23 



(2 3) (3 4) (4 5) (51) 



^scalar ^ }_^N=1 ^ ^trcc „ !££ 



[3 4] (4 1) (2 4) [4 5] ((2 3) [3 4] (4 1) + (2 4) [4 5] (5 1)) 



ln(^)-i(^-^ 

\^-S5iy 2 \s^x S23 



(3 4) (4 5) 



(3 5) [3 5f 



(12) [3 5r 



(S51 - S23,Y 

1 (12)[3 4](41)(2 4)[4 5] 



[12] [2 3] (3 4) (4 5) [51] [2 3] (3 4) (4 5) [5 1] 2 S23 (3 4) (4 5) S51 
where 



^trcc ^ ^t'■<=<=(l-^2~,3+,4+,5^ 



(12)^ 



(12) (2 3) (3 4) (4 5) (51)' 



(21) 
(22) 



20 



. ^ = r(l + 6)r^(l-e) 

^ ~ (47r)2-^ " (4^)2-^r(l - 2e) ' ^ ' 

These amplitudes contain both infrared and uhraviolet divergences, 
which have been regulated dimensionally with Z? = 4 — 2e, retaining 
terms through 0{e'^). We see that the three components have quite dif- 
ferent analytic structure, indicating that the rearrangement is a natural 
one. The N — A supersymmetric component is the simplest, followed 
by the iV = 1 component. The non-supersymmetric scalar component 
is the most complicated, yet it is still simpler than the amplitude with a 
gluon circulating in the loop, because it does not mix all three compo- 
nents together. The amplitudes for the one other helicity configuration 
needed for NLO corrections, ^5;i(l~, 2+, 3~, 4+, 5+), are a bit more 
complicated 

The surprising simplicity of TV = 4 supersymmetric loop amplitudes 
was first observed by Green, Schwarz and Brink in their calculation of 
the four-gluon amplitude as the low-energy limit of a superstring ampli- 
tude p3| . The supersymmetric decomposition can also reveal structure 
in electroweak amplitudes that would otherwise remain hidden |44[ . As 
we shall discuss in section ^, the cancellation of leading powers of loop 
momentum for supersymmetric multiplets is extremely useful for con- 
structing such amplitudes via unitarity |l4, 15 . 



4 UNITARITY 

Unitarity has been a useful tool in quantum field theory since its incep- 
tion. The Cutkosky rules |4^, ^ allow one to obtain the imaginary]^ 
(absorptive) parts of one-loop amplitudes directly from products of tree 
amplitudes. This is generally much easier than a full diagrammatic 
calculation because one can greatly simplify the tree amplitudes before 
feeding them into the calculation of the cuts. 

Having obtained the imaginary parts, one traditionally uses disper- 
sion relations to reconstruct real (dispersive) parts, up to additive ratio- 
nal function ambiguities. Although the Cutkosky rules are computation- 
ally simpler than Fcynman rules, the additive ambiguity has hampered 
their use in obtaining complete amplitudes. Here we show how this 
problem is alleviated by the supersymmetry decomposition of section ^, 
and by a complete knowledge of all functions that may enter into a 
calculation |4|, H). 

^By imaginary we mean the discontinuities across branch cuts. 
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Figure 9: The s- and t-channel cuts of a one-loop four-gluon amplitude. 
The cut lines can be gluons, fermions, or scalars. 



4.1 Cutkosky Rules 

Consider the s-channel cut of the four-point amplitude represented pic- 
torially in fig. The Mandelstam variables are as usual s = (fci -l-fe)^ 
and t — (^2 + k^)'^. According to the Cutkosky rules, the s-channel cut 
(with s > and t < 0) of this amplitude is 



iDisc A4;i(l,2,3,4) 



(27r 



\4-2f 



27r(5«(^i)27r(5W(£: 



X A*'-^^(-£i,l,2,^2)4"°( 



^2,3,4,4) 



(24) 

where ii = p and £2 = P — ki ^ k2, S^^"^ is the positive-energy branch of 
the delta-function and 'Disc' means the discontinuity across the branch 
cut. Color-ordering requires us to maintain the clockwise ordering of 
the legs in sewing the tree amplitudes. 

Suppose the amplitude had the form A^-i = cln(— s)-|-- • • = c(ln \s\ — 
in) where the coefficient c is a rational function. Then the phase 

space integral (^4|) would generate the in term but drop the In |s| term. 
Since we wish to obtain both types of terms, real and imaginary, we 
replace the phase-space integral by the cut of an unrestricted loop mo- 
mentum integral fl^ ; that is, we replace the 5- functions with Feynman 
propagators. 



A4;i(l,2,3,4) 

^4-2e 



P 



(2^)4-2. ^ A'n-ei, 1, 2, £2) ^ (-^2, 3, 4, 



(25) 
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While eq. (|2J) includes only imaginary parts, eq. ( p5| ) contains both 
real and imaginary parts. As indicated, eq. (25) is valid only for those 
terms with an s-channel branch cut; terms without an s-channel cut may 
not be correct. A very useful property of this formula is that one may 
continue to use on-shell conditions for the cut intermediate legs inside 
the tree amplitudes without affecting the result. Only terms containing 
no cut in this channel would change. A similar equation holds for the 
^-channel cut depicted in fig. Combining the two cuts into a single 
function, one obtains the full amplitude, up to possible ambiguities in 
rational functions. 



This procedure generalizes to an arbitrary number of external legs. 
Isolate the cut in a single momentum channel by taking exactly one 
of the momentum invariants to be above threshold, and the rest of 
the cyclicly adjacent ones to be negative (space-like). To construct 
all terms with cuts in an amplitude, combine the contributions from 
the various channels into a single function with the correct cuts in all 
channels. Below we describe how to link the rational functions appearing 
in amplitudes to terms with cuts, so that complete amplitudes can be 
obtained from Cutkosky rules. 



4.2 Cut Constructibility 



One-loop amplitudes satisfying a certain power-counting criterion (for 
example supersymmetric amplitudes) can be obtained directly from 
four-dimensional tree amplitudes via the Cutkosky rules. That is, when 
the criterion is satisfied, one may fix all rational functions appearing in 
the amplitudes directly from terms (through 0{e'^)) in the amplitudes 
which contain cuts. We refer to such amplitudes as 'cut-constructible'. 
(Amplitudes not satisfying the criterion can still be obtained from cuts, 
but one must evaluate the cuts to higher order in e, which is more work.) 
In the decomposition of Ar^.i given in eq. (^l|), the iV = 4 and iV = 1 
supersymmetric components are cut-constructible, while the scalar com- 
ponent is not. Correspondingly, rational functions in the first two com- 
ponents (i.e., ^ and 2) are intimately linked to the logarithms, while 
the last three rational terms in A^'if'^'-' are not so linked. 
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In a one-loop calculation one encounters integrals of the form 




(26) 



where m is the number of propagators in the loop, Ki are sums of ex- 
ternal momenta hi, and P{p^) is the loop-momentum polynomial. A 
cut-constructible amplitude is one for which one can arrange that all the 
P{p^) have degree at most m — 2, except for m = 2 when P should 
be at most linear. Any amplitude satisfying this power-counting crite- 
rion can be fully reconstructed from its cuts (through ©(e")) (l5|. The 
basic idea behind the proof is that only a restricted set of analytic func- 
tions appear in a cut-constructible amplitude. The standard Passarino- 
Veltman method |49j reduces the generic tensor integral Im[P{p^)] to 
a linear combination of basic integrals with from 2 to m external legs. 
(The kinematics of the lower-point integrals are obtained by cancelling 
denominator factors in the original integral. In a diagrammatic repre- 
sentation of the integrals, this corresponds to pinching together adjacent 
external legs.) A key feature of Passarino-Veltman reduction is that in- 
tegrals obeying the power-counting criterion can be reduced entirely to 
scalar integrals (integrals with P — 1). The proof of cut-constructibility 
is then based on showing that the cuts provide sufficient information to 
fix the coefficients of all the scalar integrals. As we shall exemplify, am- 
plitudes not satisfying the power-counting criterion contain additional 
rational functions, which spoil the argument. 

As an illustration, any cut-constructible massless four-point ampli- 
tude must be given by a linear combination of the five scalar integrals 
depicted in fig. |lO[ (The triangle integral with legs 3 and 4 pinched is 
equal to the integral with legs 1 and 2 pinched in fig. pO[ b and is there- 
fore not included in the figure; similarly, the one with 2 and 3 pinched 
is equal to the one in fig. p^.) All these integrals can be generated by 
Passarino-Veltman reduction of a box Feynman diagram; the triangle 
and bubble integrals can also be generated from other Feynman dia- 
grams. (Bubbles on external legs vanish in dimensional regularization, 
and are therefore not included.) The coefficients of the integrals are 
fixed by the cuts because each integral contains logarithms unique to 
it: the box contains the product ln(— s) ln(— t), the triangles ln(— s)^ or 
ln(— 1)2, and the two bubbles contain ln(— s) or \n{—t). Consequently 
no linear combination of these integrals with rational coefficients can be 
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w 1,1 

Figure 10: The independent scalar integrals that may appear in a 
massless four-point calculation. 



formed which is cut-free. 

The proof for an arbitrary number of external legs is similar, al- 
though more complicated. By systematically inspecting all scalar inte- 
grals that enter into an rt-point amplitude, one may show that the cuts 
fix the coefficients of all integrals uniquely [|l5) . One may also show that 
the errors induced by ignoring the difference between using D = 4 — 2e 
and D = 4 momenta on the cut legs do not affect the cuts through 
0(6°). This observation is of considerable practical use because D = 4 
tree amplitudes are simpler than those with legs in D = 4 — 2e. 

The proof breaks down for amplitudes that do not satisfy the power- 
counting criterion. For example, the scalar bubble with momentum 

miK) = J^^^^XK'X = ^ + H-K') + 2 + 0(e), (27) 

obeys the criterion. It contains a rational function, '2', but the lat- 
ter is always accompanied by ln(— ii"^). On the other hand, the linear 
combination 

{— ^-^)h[l]{KXhW]{K) = --{K'^K^-ri^'^K^)+0{e) 

(28) 

does not obey the criterion, because l2[p^p'^]{K) is quadratic in the loop 
momentum. The combination ( ^8|) is free of cuts through ©(e"); there 
is no logarithm attached to it at this order. The presence of such a 
combination within an amplitude cannot be detected using the 0(e^) 
cuts. 

In general, the power counting associated with a given amplitude 
depends on the specific gauge choice and diagrammatic organization. 
However, it suffices to find one organization of the diagrams satisfying 
the power-counting criterion. The string-inspired method discussed in 
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section|^ provides such an organization; it can satisfy the power-counting 
criterion even when the corresponding diagrams in conventional Feyn- 
man gauge do not. An important class of cut-constructible amplitudes 
are those in supersymmetric gauge theory. In section 3.2 we showed that 
for n-gluon amplitudes the leading two powers of loop momentum cancel 
in a supermultiplet contribution; the same result holds for amplitudes 
with external fcrmions . 



4.3 Supersymmetric Examples 

As a simple example, consider the contribution of an = 4 super- 
symmetry multiplet to a four-gluon amplitude. This amplitude is an 
ordinary gauge-theory amplitude but with a particular matter content: 
one gluon, four gluinos and six real scalars all in the adjoint representa- 
tion. As discussed in section ^ 2+, 3+, 4+) = so the first 
non-trivial case to consider is A^.^'*'{1~ , 2~ , S"*", 4+). 

For the s-channel cut depicted in fig. ^, only the gluon loop con- 
tributes; for ferniion or scalar loops the supersymmetry identities in 
eq. (|lj) guarantee that at least one of the two tree amplitudes vanish. 
The necessary tree amplitudes are the four-gluon amplitudes 



(12^"^ 

Al"'%-i+,i-,2-,e+) = 1- ^ ' 



itsaci n- o+ j+ : \ -^1 



v4 



(29) 



-f2 3)(3 4)(4£i) {l^ ~i2 



All other combinations of helicities of the intermediate lines cause at 
least one of the tree amplitudes on either side of the cut to vanish. (The 
outgoing-particle helicity convention means that the helicity label for 
each intermediate line flips when crossing the cut.) Cut-constructibility 
of supersymmetric amplitudes allows us to use the four-dimensional tree 
amplitudes, so that the cut in the s channel, eq. (pq), becomes 



<r'(l-,2-,3+,4+) = / 

s-cut J 



(27r)4-2e ^2 (£,l)(12)(2^2)(^2^l) 



il (^2 3)(3 4)(4^i)(4^2 

(30) 

where we have removed the minus signs from inside the spinor products 
by cancelling constant phases. To put this integral into a form more 



s— cut 
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reminiscent of integrals encountered in Feynman diagram calculations 
we may rationalize the denominators using, for example, 



[2^2 



(2 £2) ip-h)^ 



(31) 



We use the on-shell conditions — and £2 = 0, which apply even 
though the loop integral is unrestricted, because of the s-cut restriction. 
Performing such simplifications yields, 



<r'(l",2-,3+,4+) 



S-cut 



(27r)4-2^ p2(p _ fc^)4(p _ fc^ _ ^2)2(p + k^Y 

where we have extracted a factor of the tree amplitude, 



Ar''(l-,2-,3+,4+) = z- 



(12) 



(12) (2 3) (3 4) (41) 
from the amplitude. The numerator of the integrand is 

M = [h 1] (1 4) [4 (4 12) [4 3] (3 2) [2 £2] (^2 h) 
= tr+[£il4£i432£2] 

= -4tr+[4321] ^1 • fci ^1 • fc4 = -St (p - kif{p + k^f , 



S-cut 

(32) 



(33) 



(34) 



where tr+[- • •] = ^tr[(l + 75) • • •] and we used 



0, 



+ (35) 



The 75 term in the trace drops out because a four-point amplitude 
has only three independent momenta to contract into the totally anti- 
symmetric Levi-Civita tensor. 

Thus in eq. ( ^2|) the numerator neatly reduces the squared propaga- 
tors to single propagators, 



d 



4-2, 



P 



1 



(36) 



(2^)4-2e p2(p _ fc^)2(p _ fc^ _ fc2)2(p + k^f 

which is a scalar box integral. Thus the s-cut contribution is given by 

-st 



Air(l-,2-,3+,4+) 



s-cut (47r)2 



(37) 
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where the massless scalar box integral is (see e.g. ref. psf) 
/4(.,i) = -^|-^[(-.)-^ + H)-^]+iln2(^)+y|+0(6). (38) 

The evaluation of the i-channel cut depicted in fig. ^ is similar, but 
a bit more involved since all particles in the multiplct contribute. How- 
ever, after summing over the contribution of all particles, with the help 
of the SWI (^4|) and the Schouten identity (|Tl|), the integral appearing 
in the f-channel cut turns out to be the same as the one appearing in 
the s-channel cut in eq. (^). 

Combining the s and t channel results, the amplitude must be 

A^r'(l",2-,3+,4+) ^ ^—^4-/4(5, t). (39) 

The rational function proportional to tt^ contained in the box integral 
( |3^ ) is fixed by the cuts since it appears in association with the loga- 
rithms in this function. Integrals having cuts in multiple channels, such 
as 74(5, i), provide a strong consistency check: their coefficients can be 
obtained via two or more separate cut calculations and the results must 
agree. 

Following the same procedure one may evaluate the other nonvanish- 
ing N = A four-gluon amplitude, A^.^'^{1~ , 2+, 3~, 4+), where the nega- 
tive helicities are non-adjacent. Surprisingly, the same basic calculation 
can be easily extended to an arbitrary number of external legs for max- 
imally helicity violating (MHV) amplitudes, those with two negative- 
helicity gluons and the remaining of positive helicity. (A special case 
is ^^f4(l-,2-,3+,4+,5+), given in eq. (|l|).) The cuts faU into two 
categories, depending on whether the external negative-helicity gluons 
are on the same or on opposite sides of the cut, as depicted in fig. |l^. 
In cither case the tree amplitudes on both sides of the cuts are given by 
the Parke- Taylor formula pu, H] , 



A^'^'^'^ {l^ , nil , . . . , k '■■•jJ T ■ ■ Tf^2 T^'t) 



{£1 mi) (mi, mi -1-1) • • • (m2 — 1, TO2) (m2^2) (^2^1) 

(40) 

where j and k are the two negative-helicity legs, or by formulae related 
to eq. (|^ by the SWI (^. The key to evaluating the cut integrals 
for an arbitrary number of external legs is that only two denominator 
factors in the tree amplitudes ( ^0[ ) contain the loop momentum (since 
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Figure 11: The relevant cuts for computing the MHV amphtudes for 
an arbitrary number of external legs. 




Figure 12: AU-n MHV supersymmetric amplitudes can be evaluated 
by evaluating hexagon integrals. 

1/ {(-2 ^i) — [^2 ^i] / (fcmi + • • • + km^Y). Thus each tree contributes only 
two propagators containing the loop momentum, so after including the 
two cut propagators the hardest integral to be evaluated is the hexagon 
integral depicted in fig. |lj. These hexagon integrals can be reduced to 
scalar box integrals in much the same way as for the four-point case, al- 
lowing one to obtain the amplitudes for an arbitrary number of external 
legs 0. 

The analysis of = 1 supersymmetric MHV amplitudes is similar, 
although more complicated |]l5j. Again the key to the construction is 
that no more than six denominators contain loop momentum, even for 
an arbitrary number of external legs. One instance of the general N = 1 
MHV result is provided by A^f 1(1", 2", 3+, 4+, 5+) in eq. (|l|). Notice 
that only the S23 and S51 channels contain cuts. This result (which is 
also true for the scalar component) is a simple consequence of the super- 
symmetry identities ( p^ ) . The construction of amplitudes via cuts does 
not rely on supersymmetry, but only on the power-counting criterion; 
however, non-supersymmetric amplitudes generally do not satisfy the 
criterion. 
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4.4 N on- super symmetric Example 

Amplitudes not satisfying the power-counting criterion require an ex- 
tension of this approach. Consider the non-supersymmetric amphtude 
for four identical helicity gluons with a scalar in the loop, 

A|--(1^2+,3+,4+) = -^^^m (41) 

first obtained from string-based techniques |Q . At first sight one might 
think that it is impossible to use unitarity to obtain this amplitude, 
since it contains no cuts. The box Feynman diagram for this amplitude 
contains up to four powers of loop momentum, so the power-counting 
criterion is not satisfied (in any gauge). 

However, in D = 4 — 2e all terms in a massless amplitude necessar- 
ily have cuts by dimensional analysis of eq. (^6|), all terms must 
be proportional to factors of (— X^)"'^, where is some kinematic 
variable. In particular a massless four-point amplitude must be of the 
form 

= (1 - eln(-s))/i + (1 - eln(-t))/2 + • • • 

where /i and /2 are dimensionless functions of the kinematic variables. 
This expression now contains cuts at 0{e) even if /i and /2 are cut-free. 
Rational functions such as those in eq. may therefore be obtained 
from the sum /i + /2, fixed by the coefficients of the single logarithms 
at 0(e). 

Thus, to obtain the rational function contributions in amplitudes 
which do not satisfy the power-counting criterion we must perform a 
cut calculation valid to at least one higher order in e. We are not 
actually interested in the explicit values of the 0(e) terms; we only 
need to extract the sum fi + f2- To implement a calculation valid 
to higher orders in e we correct for the fact that the loop momenta 
appearing in the tree amplitudes on either side of the cut are in (4 — 2e)- 
dimensions instead of four-dimensions. The proper on-shell conditions 
on the cut legs are £1 — — and ^2 ~ = 0, where £1 and £2 
are left in four-dimensions and fi is the (— 2e)-dimensional part of the 
loop momentum. We follow the standard prescription that the (— 2e)- 
dimensional subspace is orthogonal to the four-dimensional one [g9| . For 
practical purposes we may think of as a mass which gets integrated 
over. (This decomposition of the loop momentum has also been used 
by Mahlon in his recursive approach.) 
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For the amplitude in eq. (41), the tree amplitudes entering the two 
sides of the s-channel cut, depicted in fig. are easily computed from 
color-ordered Feynman diagrams; the one on the left side of the cut is 



where legs £i and £2 represent the cut scalar lines. The one on the 
right side is obtained by relabeling legs. The amplitude ( ^3|) vanishes in 
D — A {y? "> 0) by the SWI (p^). Plugging these tree amplitudes into 
eq. (pq), we obtain the s cut of the scalar loop contribution. 



^scalar(^+^2+,3+,4^ 



where the factor of 2 is from the two states of a complex scalar and 



(27r)4 (27r)-2- (p^ - ^i^){{p-k^Y ^ f^'') ■ ■ ■ {{p+U? - A^^) ' 

(45) 

The t-channel cut, depicted in fig. is similar and may be obtained 
via the relabeling 1 <-> 3. Using the identity 

[3 2] [14] _ [12] [3 4] ^^g^ 



(3 2) (14) (12) (3 4)' 

the i-cut is given simply by eq. (p^, with 's-cut' replaced by 't-cut'. 
Combining the two cuts we obtain an expression valid for both cuts, 

Alf-il+XXA-") = 2|^|||^4[/.1 . (47) 

Although we only calculated the cuts, we did so to all orders in e; there- 
fore by eq. ( ^2| ) we know the complete loop amplitude. To obtain the 
amplitude through 0{e^) we need only evaluate the leading ©(e") con- 
tribution to the integral X4 [/i'*] . 

A good way to evaluate the leading term is to first integrate out 
the angles in the (— 2e)-dimensional subspace. Using the fact that the 
integrand is a function only of /i^, we have 

^ ■ -(4-)^7.7A-T / V(M^)-^-^ (48) 



(27r)-2<^ ' ' r(l 

The overall e from the measure must be compensated by a 1/e ultraviolet 
pole in the remaining integration. As usual a leading ultraviolet diver- 
gence may be extracted conveniently by setting all external momenta to 
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zero and inserting a mass parameter A to replace the momenta. Thus 
we may evaluate the integral ( ^5| ) as 



7 (27r)4 (2^)-2^ (p2_^2_^2)4 

where we have used standard formulas |46| for the angular integrals and 
then integrated the radial dimension. Plugging the leading-in-e result 
into eq. (^), we obtain the correct result for the amplitude (^). 

Although this method can in principle be applied to any massless 
one-loop amplitude to obtain complete amplitudes, it is generally ad- 
vantageous to first decompose amplitudes into pieces which are cut- 
constructible and pieces which are not. One may also calculate loop 
amplitudes for massive particles in this way [Q, but cut-free integrals 
may appear. The coefficients of these functions must be determined by 
other means, such as knowledge of ultraviolet or infrared divergences. 



5 FACTORIZATION 

In quantum field theory, amplitudes are constrained by their behavior 
as kinematic variables vanish; they must factorize into a product of two 
amplitudes with an intermediate propagator. This may be used as a 
check on five- or higher-point amplitudes. (Factorization of four-point 
amplitudes in a theory without massive particles is trivial since the 
limiting kinematics is degenerate.) Factorization properties may also be 
used to help construct new amplitudes from known ones. In principle, 
this can be an extremely efficient way to obtain amplitudes since one 
avoids evaluating loop integrals. 

Mangano and Parke have reviewed the factorization properties of 
tree-level QCD amplitudes ||]. We shall focus on the corresponding 
properties at one loop, which are a bit more complicated since the am- 
plitudes generally contain infrared divergent pieces which do not factor- 
ize naively. Nevertheless, as any kinematic variable vanishes, one-loop 
amplitudes have a universal behavior quite similar to that of tree-level 
amplitudes. 



32 



5.1 General Framework 



First we review briefly the situation at tree level. Color-ordered am- 
plitudes can have poles only in channels corresponding to the sum of 
cyclicly adjacent momenta, that is as P^j — > 0, where P^j = (ki + 
ki+i + • • • + fcj)^. This is because singularities arise from propagators 
going on-shell, and propagators for color-ordered graphs always carry 
momenta of the form P^^ . The general form of an ri-point color-ordered 
tree amplitude in the limit that Pi^ vanishes is 

^ 7 (50) 

^ Al^^l, . . . , m, P>^) AT.'^_^M -I- 1, . . . , n, p-^) , 

where Pi.m is the intermediate momentum, A*-^^^ and A*-^^^_^t^ are 
lower-point scattering amplitudes, and A denotes the helicity of the in- 
termediate state P. The intermediate helicity is reversed going from one 
product amplitude to the other because of the outgoing-particle helicity 
convention. 

For two-particle channels (m = 2), eq. ( |50| ) needs to be modified, 
because a three-point massless scattering amplitude is not kinematically 
possible. As PI2 0, ki and ^2 become collinear. QCD amplitudes have 
an angular-momentum obstruction in this limit. For example, a gluon 
of helicity +1 cannot split into two collinear helicity ±1 gluons and 
conserve angular momentum. This transforms the full pole in P^2 = S12 
into the square-root of a pole, 1 / y'si2, a behavior which is well captured 
via the spinor products (12) ,[12]. It is useful to lump all terms not 
associated with A^jl2.\ in eq. ( pO| ) into a 'splitting amplitude' Split*''''''. 
In particular, as the momenta of adjacent legs a and b become collinear, 
we have 

<-(. . . , a^", . . .) ^ Split*-T(2, ^'^^) ■■,P\---), 

(51) 

where P is the intermediate state with momentum kp = ka + kb, A 
denotes the helicity of P, and z is the longitudinal momentum fraction, 
ka ~ zkp, fcf, « (1 — z)fcp. The universality of these limits can be derived 
diagrammatically, but an elegant way to derive it is from string theory 

because all the field theory diagrams on each side of the pole are 
lumped into one string diagram. 

Given the general form ( pT] ) , one may obtain explicit expressions for 
the tree-level g gg splitting amplitudes from the four- and five-gluon 
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amplitudes ^, |2l| . For example, taking the coUinear limits of eq. 
for ^^'■'"'(l^, 2^, 3+, 4+, 5+) and comparing to eqs. ( |33| ) and (51) shows 
that 



Split*!^"(a-,6-) = 0, 
Split*!™ (a+, 6+) 

Split*;°'=(a+,6-) 

Split*;°'=(a-,fe+) 



^/z{l-z) {ab) ' 
{ab) ' 



^z{l-z) (ab) 



(52) 



The remaining helicity configurations are obtained using parity. The 
g qq and q qg splitting amplitudes can be obtained in similar 
fashion. 

The situation for color-ordered one-loop amplitudes is similar to tree 
level. The one- loop analog of eq. (|5^) is schematically depicted in fig. |l^, 
and is given by 



l,m 



A'^^;,{1, ...,m, P^)'A':Z.+iim + l,...,n, P"^) 



(53) 

where the one- loop factorization function, Fact„, is independent of he- 
licities and does not cancel the pole in In an infrared divergent 

theory, such as QCD, amplitudes do not factorize 'naively': Factn may 
contain logarithms of kinematic invariants built out of momenta from 
both sides of the pole in P^^^; ln(— s„4) is an example of such a loga- 
rithm. The factorization functions are nonetheless universal functions 
depending on the infrared divergences present in the amplitudes p4[ |. 

The collinear limits for color-ordered one-loop amplitudes are a spe- 
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Figure 13: A schematic representation of the behavior of one- loop 
amplitudes as a kinematic invariant vanishes. 






b 

Figure 14: A schematic representation of the behavior of one-loop 
amplitudes as the momenta of two legs become collinear. 

cial case and have the form 

^loop ^ Y.fspl\t'lTiz,a'-,b^^)A':°li...ia + b)\..) 

(54) 

+ SplitLT(^, a^\b^') ATii- . . (a + . .) ) , 



which is schematically depicted in fig. 14, The splitting amplitudes 
Spht!!^°(a'^° , b'^'') and Split'f^''(a'''" , b^'') are universal: they depend only 
on the two momenta becoming collinear, and not upon the specific am- 
plitude under consideration. The explicit Split'°^''(a^" , 6'*''') were orig- 
inally determined from the four- and five-point one-loop amplitudes 
p^ , in much the same way as we obtained the tree-level splitting 
amplitudes above. (See appendix B of ref. Soft limits — the be- 

havior as any particular fc^ — > — are also useful for constraining the 
form of one-loop amplitudes, and have a form analogous to eq. (p^. 

In performing explicit calculations, factorization provides an ex- 
tremely stringent check since one must obtain the correct limits in all 
channels. A sign or labeling error, for example, will invariably be de- 
tected in some limits. In some cases one can also use factorization to 
construct ansatze for higher-point amplitudes Q . One writes down 
a sufficiently general form for a higher-point amplitude, containing arbi- 
trary coefficients which are then fixed by imposing the correct behavior 
as kinematic variables vanish. A collinear bootstrap of this form would, 
however, miss functions that are nonsingular in all collinear limits. For 
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five-point amplitudes it is possible to write down such a function, namely 



(12) (2 3) (3 4) (4 5) (51) ' ^ ' 

since the contracted antisymmetric tensor e(l, 2, 3, 4) = ^iSfj^up^k'^kl^k^kl 
vanishes when any two of the five vectors ki become collinear (using 
X]i=i ki — Q). However, it is quite possible that the factorization con- 
straint uniquely specifies the rational functions of color-ordered n > 6- 
point amplitudes, given the lower point amplitudes. A heuristic expla- 
nation of this conjecture is that as the number of external legs increases, 
by dimensional analysis the amplitudes require ever increasing powers of 
momenta in the denominators. Thus one expects more kinematic poles 
from the denominator than zeros from the numerator. We know of no 
counter-examples to this conjecture, but don't have a proof either. 

5.2 Examples 

As an example of the behavior of a one-loop amplitude in a collinear 
limit, consider the = 4 five-gluon amplitude A^.^^ given in eq. (|2l|). 
Taking the limit fc4 — > zP, k^, ^ (1 — z)P, and using the four-gluon 
result fp9h, we find 



A^r^(l-, 2-, 3+, 4+, 5+) ^ Split*!^^(4+, 5+)A^r^(l-, 2-, 3+, P+) 

+ Split^=4(4+, 5+X°°(l-, 2-, 3+, P+) , 

_ (56) 

where the tree splitting amplitude is given in eq. (p3). This limit deter- 
mines the one-loop A^ = 4 multiplet contribution to the g gg splitting 
amplitude. 



Split^=^(a+,6+) 



cr Split";''(a+,6+) 



^-{-Sabz{l - z))-'=-H21nzln(l - z)- — 



1 

V- -/ g 

(57) 

Note that the loop splitting amplitude has absorbed the mismatch of 
infrared divergences between the five- and four-point amplitudes. This 
splitting amplitude will reappear in different collinear limits of this and 
other amplitudes. 

Given the splitting amplitudes and five-point amplitudes, it is also 
possible to construct conjectures for higher-point amplitudes by de- 
manding that they factorize correctly. Consider, for example, the complex- 
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scalar loop contribution to a five-gluon amplitude with all identical he- 
licities |||, 

^scalar(^+^2+,3+,4+,5 + ) 

_ i 512S23 + 523^34 + ^34^45 + ^45^51 + 551^12 + g(l, 2, 3, 4) 

~ 967r2 (12) (2 3) (3 4) (4 5) (51) 

(58) 



As noted in section 3.2, for this helicity configuration the gluon and 
fermion loops are proportional to the scalar-loop contribution. One can 
verify that this amplitude has the correct coUinear limits (p^, using the 
four-gluon amplitude (^). 

Using eq. (^J), the explicit form of the tree splitting amplitudes 
(|5^, yl'[°°(l*, 2+, • • • , n+) = 0, and experimenting at small n, we can 
construct higher-point amplitudes by writing down general forms with 
only two-particle poles, and requiring that they have the correct coUinear 
limits. Doing so leads to the all-n ansatz [ p^ , 

.scalar/. + ^+ „ + N ^ ^_ -[^2^^] 

^n;l 1^ , ^ ) 43^3 (1 2) (2 3) • • • (u 1) ' 

(59) 

where tr_ [iii2i3j4] = ^tr_[(l — ^^)^i^^j^^'^i^^iA.Th\s result has been 
confirmed by Mahlon via recursive techniques p3[ . 

Indeed, the infinite sequence of one-loop iV = 4 supersymmetric 
MHV amplitudes was first constructed via a coUinear bootstrap and 



only then calculated using the unitarity method described in section|4.3. 
Other helicity configurations are more complicated, due to the appear- 
ance of multi-particle poles. Nevertheless, one can construct some six- 
point amplitudes from knowledge of the five-point amplitudes. This is 
most useful for the rational-function parts, which can be obtained via 
unitarity only by working to higher order in e. 



6 CONCLUSIONS AND OUTLOOK 

We have reviewed various developments in calculational techniques for 
one-loop gauge theory amplitudes, especially in QCD. Such calculations 
are necessary in order to confront theoretical predictions with experi- 
ments to some degree of precision. Feynman rules, however, become 
extremely cumbersome for one-loop multi-parton calculations. Even the 
simplest processes are rather difficult to calculate without aid of a com- 
puter and for five or more external legs traditional methods break down 



37 



because of an exponential explosion in algebra. The results, however, 
are usually quite compact, especially when compared to intermediate 
expressions. 

The computational situation can be greatly improved by combining 
a number of ideas. Methods that have previously been used at tree-level, 
such as spinor helicit y , color decomposition , and supersymmetry 
Ward identities [^8|, |19[, remain very useful at one loop. String theory 
motivates a number of improved organizational ideas such as supersym- 
metric decompositions, relations between color-decomposed amplitudes 
and improved gauge choices |3^ , |3^ , These ideas mesh nicely 
with the use of Cutkosky rules iHSj to obtain complete amplitudes. In 
the superstring organization of the amplitude, components can be iden- 
tified whose rational as well as cut-containing parts can be obtained 
directly from knowledge of the branch cuts The remaining 

components, though more difficult, can be attacked either by evaluat- 
ing cuts to higher order in e or by exploiting universal factorization 
properties |p^ . 

The techniques discussed in this review have made possible a va- 
riety of new calculations, including those of all five-parton amplitudes 
, |2^, Q and of certain infinite sequences of massless amplitudes. The 
methods have also been applied to amplitudes containing massive parti- 
cles m and to gravitational amplitudes Q . Mahlon has also used 
recursion relations, outside of the scope of this review, to obtain infinite 
sequences of fermion loop amplitudes with maximal helicity violation 

It would be desirable to extend these techniques to two-loop multi- 
parton calculations; while various authors have taken first steps |5^, 
|55| , |5^ in this direction, a great deal of work remains to be done. 

We thank G. Chalmers, A.G. Morgan and especially D.C. Dunbar 
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ported in part by the US Department of Energy under grants DE-FG03- 
91ER40662 and DE-AC03-76SF00515, by the Alfred P. Sloan Foun- 
dation under grant BR-3222, and by the Direction des Sciences de la 
Matiere of the Commissariat d I'Energie Atomique of France. 
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